Light WIMP search in XMASS 



K. Abe a,c , K. Hieda a , K. Hiraide a,c , S. Hirano a , Y. Kishimoto a,c , 
K. Kobayashi a ' c , S. Moriyama a,c , K. Nakagawa a , M. Nakahata a ' c , 
H. Ogawa a,c , N. Oka a , H. Sekiya a,c , A. Shinozaki a , Y. Suzuki a ' c , 
A. Takeda a ' c , O. Takachio a , K. Ueshima a ' 1 , D. Umemoto a , M. Yamashita a,c , 
B. S. Yang a , S. Tasaka b , J. Liu c , K. Martens c , K. Hosokawa 6 , K. Miuchi 6 , 
A. Murata 6 , Y. Onishi 6 , Y. Otsuka 6 , Y. Takeuchi e ' c , Y. H. Kim f , K. B. Lee f , 
M. K. Lee f , J. S. Lee f , Y. Fukuda s , Y. Itow h ' d , K. Masuda h , Y. Nishitani h , 
H. Takiya h , H. Uchida h , N. Y. Kim 1 , Y. D. Kim 1 , F. Kusaba j , D. Motoki^ 1 , 
K. Nishijima-', K. Fujii 1 , I. Murayama 1 , S. Nakamura 1 

a Kamioka Observatory, Institute for Cosmic Ray Research, the University of Tokyo, 
Higashi-Mozumi, Kamioka, Hida, Gifu, 506-1205, Japan 
b Information and multimedia center, Gifu University, Gifu 501-1193, Japan 
c Kavli Institute for the Physics and Matematics of the Universe, the University of 

Tokyo, Kashiwa, Chiba, 277-8582, Japan 
d Kobayashi-Maskawa Institute for the Origin of Particles and the Universe, Nagoya 
University, Furo-cho, Chikusa-ku, Nagoya, Aichi, 464-8602, Japan. 
e Department of Physics, Kobe University, Kobe, Hyogo 657-8501, Japan 
f Korea Research Institute of Standards and Science, Daejeon 305-340, South Korea 
9 'Department of Physics, Miyagi University of Education, Sendai, Miyagi 980-0845, 

Japan 

h Solar Terrestrial Environment Laboratory, Nagoya University, Nagoya, Aichi 464-8602, 

Japan 

l Department of Physics, Sejong University, Seoul 143-747, South Korea 
J Department of Physics, Tokai University, Hiratsuka, Kanagawa 259-1292, Japan 
k School of Science and Technology, Tokai University, Hiratsuka, Kanagawa 259-1292, 

Japan 

l Department of Physics, Faculty of Engineering, Yokohama National University, 
Yokohama, Kanagawa 240-8501, Japan 



Abstract 

A search for light dark matter using low-threshold data from the single 
phase liquid xenon scintillation detector XMASS, has been conducted. Using 
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the entire 835 kg inner volume as target, the analysis threshold can be lowered 
to 0.3keVee (electron-equivalent) to search for light dark matter. With low- 
threshold data corresponding to a 5591.4 kg-day exposure of the detector and 
without discriminating between nuclear-recoil and electronic events, XMASS 
excludes part of the parameter space favored by other experiments. 
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1. Introduction 

There is substantial evidence that the universe is composed of dark mat- 
ter [l[ . Among the most plausible dark matter candidates, Weakly Interact- 
ing Massive Particles (WIMPs), which are expected to couple to ordinary 
matter primarily through the weak force, can be detected directly through 
observation of nuclear recoils produced in their elastic scattering interac- 
tions with detector nuclei (2)]. Many theories of physics beyond the Standard 
Model predict stable WIMPs with supersymmetric models still among the 
most popular In the constrained Minimal Supersymmetric Standard 
Model, the lightest supersymmetric particle is favored as a WIMP candidate 
and has an expected mass of M x > 100 GeV. However, some experiments 
indicate a possible WIMP signal 0] with a lighter mass of ~ 10 GeV 
and with a spin- independent cross section of the order of ~ 10 _40 cm 2 . These 
positive signals have come predominantly from experiments without the ca- 
pability to discriminate between electromagnetic and nuclear recoils, while 
other experiments that have this ability have excluded light WIMPs at these 
cross sections 0, 0, 0] ■ This Letter presents a search for light WIMPs using 
a 5591.4 kg-day exposure of the XMASS experiment without nuclear recoil 
discrimination. 



2. XMASS detector 

The XMASS experiment is located underground in the Kamioka Obser- 
vatory at a depth of 2700 m.w.e. and has been designed to search for WIMP 
interactions using the self-shielding properties of liquid xenon in an ultra low 



radioactivity environment 10]. XMASS is a single phase liquid xenon scin- 
tillator detector containing 1050 kg of Xe in an OFHC copper vessel. Xenon 
scintillation light is collected by 630 hexagonal and 12 cylindrical inward- 
pointing Hamamatsu R10789 series photomultiplier tubes (PMTs) arranged 
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on an 80 cm diameter pentakis-dodecahedron support structure within the 
vessel to give a total photocathode coverage of 62.4%. 

These PMTs view an active target region containing 835 kg of liquid 
xenon. In order to monitor the PMT stability and measure the trigger ef- 
ficiency, eight blue LEDs with Teflon diffusers are mounted to the support 
structure. There are six LEDs arranged along the equator and one each at the 
top and the bottom of the pentakis-dodecahedron. To shield the scintillator 
volume from external gammas, neutrons, and muon-induced backgrounds, 
the copper vessel is placed at the center of a 010 m x 11m cylindrical tank 
filled with pure water. This volume is viewed by 72 Hamamatsu R3600 20- 
inch PMTs to provide both an active muon veto and passive shielding against 
these backgrounds. This is the first water Cherenkov shield used in a dark 
matter search experiment. To perform energy and position reconstruction 
calibrations a portal has been prepared along the central vertical axis (z- 
axis) of the PMT support structure through which an OFHC copper rod 
can be inserted directly into the target volume. Thin cylindrical calibration 
sources containing either of 55 Fe, 57 Co, 109 Cd, or 241 Am are placed at the tip 
of this rod to perform detector calibrations. A more detailed description of 
the XMASS detector is presented in [ill ]. 

PMT signals are passed though preamplifiers with a voltage gain factor 
of 11 before being processed by Analog-Timing-Modules (ATMs) [12j]. These 
modules combine the functions of typical ADC and TDC modules, recording 
both the integrated charge and the arrival time of each PMT signal. For each 
PMT channel the discriminator threshold is set to — 5mV, which corresponds 
to 0.2 photoelectrons (p.e.). When a PMT signal is above threshold a "hit" is 
registered on the module. A global trigger is generated based on the number 
of hit PMTs within a 200 ns window. 

A complete XMASS detector Monte Carlo (MC) simulation package based 



on Geant4 [_13y and including readout electronics has been developed [11 
and used for this analysis. The simulation has been tuned using calibration 
data and the optical properties of the liquid xenon have been determined 
empirically using data taken at nine points along the z-axis: z = ±40 cm, 
±30 cm, ±20 cm, ±10 cm, cm. 

3. Data and Event Selection 

The data used for this analysis, corresponding to 6.70 days of livetime, 
was taken in February 2012 with a low trigger threshold of four PMT hits. 
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Using 122 keV gammas from the 57 Co calibration source the xenon light yield 
was found to be 14.7photoelectrons/keVee. This large light yield allows 
the analysis threshold to be lowered sufficiently for sensitivity to low mass 
WIMPs. In order to achieve optimal sensitivity, the entire detector volume 
is used because fiducialization is increasingly difficult at these low energies. 

In the analysis presented below, the nuclear-recoil equivalent energy (E NR ) 
is determined by conversion from electron equivalent energy (E ee ) using the 
scintillation efficiency, £ e g, for nuclear recoils relative to that of 122 keV 
gammas at zero electric field. The treatment of the energy dependence of 



CeS and its uncertainty is the same as in Ref. [14|. The resulting functional 
dependence of E NR on E ee with its one-cx uncertainty is shown in Fig. [TJ To 
illustrate the impact of this uncertainty on the present analysis, the 0.3 keVee 
analysis threshold is also shown in the figure. In the £ e ff framework, though 
the absolute energy scale is determined at 122 keVee, dark matter signals are 
expected to appear at lower energies. Lower energy calibration samples are 
used to evaluate the accuracy of the detector simulation at these lower ener- 
gies. Fig. [2] shows the ratio of the total number of observed photoelectrons 
relative to the simulation's prediction as a function of the electron equivalent 
energy for several calibration samples. The apparent deviation of the sim- 
ulation from the data reflects not only imperfection in the modeling of the 
detector environment but also the non-linear response of electron-mediated 
events in liquid xenon [15]. Since the 5.9 keV X-ray from 55 Fe is the lowest 
energy calibration point, the response at lower energies is extrapolated using 
a linear fit through all calibration energies. In the analysis the fitting error 
from this procedure, 1.3% at 0.3keVee, is treated as a systematic effect on 
the energy scale. However, the effect of this uncertainty is small relative to 
that induced by the uncertainty in £ eS , which is 13 % at the same energy. 

The trigger efficiency as a function of the total number of photoelec- 
trons near the trigger threshold is smeared by the response of the electron- 
ics. Moreover, the trigger efficiency depends upon the WIMP mass since 
the distribution of the expected number of hits is a function of that mass. 
For this reason the trigger efficiency near threshold is derived using LED 
data taken at single photoelectron level intensities in conjunction with the 
detector simulation. During the measurement LEDs are flashed at 100 Hz 
and each flash of the LED is recorded regardless of the number of PMT hits 
observed. Based on this information the distribution of the number of hit 
PMTs, the nhit distribution, is formed for each LED. Using the detector 
simulation the same distribution is made for each WIMP mass (The astro- 
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physical parameters used in this simulation are described in Section 2]). Each 
LED event is then reweighted based on the ratio of the nhit distributions so 
that the number of expected photoelectrons, and hence energy, for a WIMP 
signal can be derived from the LED data. Comparison of the number of LED 
flashes with the number of global triggers issued by the ATMs is then used to 
compute the expected WIMP trigger efficiency as a function of the observed 
energy. The resulting efficiency curves for representative masses from the 
allowed region of ji[ are shown in Fig. |3j The analysis threshold is chosen 
as the energy at which the trigger efficiency is greater than 50 % for 5 GeV 
WIMPs and corresponds to 0.3keVee. The error bars in the figure illustrate 
that the trigger efficiency varies somewhat based upon which LED is used. 
This variation reflects intrinsic differences in the LEDs themselves as well 
as a possible position dependence of the detector response and is therefore 
treated as a systematic error in the analysis. A 1.5% error is determined 
at 0.3keVee. Using only the predicted nhit distribution from the WIMP 
simulation, the efficiency curve's shape is well reproduced within errors. 

The energy spectrum for all events in the entire volume of the detector 
over a 5591.4 kg-day exposure is shown in Fig. [5J To remove backgrounds 
prior to analysis three types of timing cuts are applied to the raw data. 
In order to remove events caused by the tail of the scintillation light dis- 
tribution observed in energetic gamma, beta, alpha, and muon interactions 
within XMASS, events that occurred within 10 ms of the previous event are 
rejected. Events whose timing distribution has an RMS greater than 100 ns 
are also removed. The last and most significant cut is designed to remove 
events which produce light predominantly through Cherenkov emission, in 
particular the beta particle emitted in the decay of 40 K contamination in 
the PMT photocathode. To select these events a time-of-flight correction is 
made to the timing distribution of each event assuming the event vertex is 
on the surface of the PMT closest to the charge-weighted center of gravity 
of the event. After this correction the timing distribution of Cherenkov-like 
events is found to be narrower than that for scintillation-like events. Events 
with more than 60 % of their PMT hits occurring within the first 20 ns of 
the event window are removed as Cherenkov-like. The ratio of the number 
of PMT hits within the first 20 ns relative to the total number of hits in the 
event window for all events (head-to-total ratio) is shown in Fig. HJ Each 
step of the data reduction is shown in Fig. |5j 

The expected WIMP acceptance efficiency of these cuts was estimated 
with the detector simulation. In the simulation WIMP recoil energy spec- 
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tra were generated for each WIMP mass and MC events were distributed 
uniformly throughout the detector volume using a liquid scintillation decay 
constant of 25 ns (l6| . Fig. [6] shows the resulting signal acceptance efficiency 
at energies below 1 keVee. The size of the error bars comes primarily from the 
systematic uncertainty in the xenon scintillation decay constant, 25 ± Ins, 
which is estimated based on the difference between the XMASS model 16 
and the NEST model 17| based on 18] . A systematic error on the selection 



efficiency is determined based on the error resulting from a linear fit to the 
points in the figure. At the 0.3keVee analysis threshold this error is 6.1%. 



4. Results and Discussion 

Figure [7] shows simulated WIMPs energy spectra overlaid on the observed 
spectrum after the data reduction was applied. WIMPs are assumed to 
be distributed in an isothermal halo with v Q = 220km/s, a galactic escape 
velocity of t> esc = 650km/s, and an average density of 0.3 GeV/cm 3 . In order 
to set a conservative upper bound on the spin-independent WIMP-nucleon 
cross section, the cross section is adjusted until the expected event rate in 
XMASS does not exceed the observed one in any energy bin above 0.3keVee. 
Implementing the systematic errors discussed in the text above, the resulting 
90 % confidence level (C.L.) limit derived from this procedure is shown in 
Fig. [HI The impact of the uncertainty from £ e ff is large in this analysis, so 
its effect on the limit is shown separately in the figure. 

After careful study of the events surviving the analysis cuts, their origins 
are not completely understood. Contamination of 14 C in the GORE-TEX® 
sheets between the PMTs and the support structure may explain a fraction 
of the events. Light leaks through this material are also suspect. Nonethe- 
less, the possible existence of a WIMP signal hidden under these and other 
backgrounds cannot be excluded. Although no discrimination has been made 
between nuclear-recoil and electronic events, and many events remain in the 
analysis sample, the present result excludes part of the parameter space fa- 
vored by other measurements [1, 0, 0] when those data are interpreted as a 
signal for light mass WIMPs. Finally, we are working on modifications to 
the inner surface of XMASS, especially around the PMTs, to improve the 
detector performance. 
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Figure 1: The relationship between the electron equivalent energy and the nuclear-recoil 
equivalent energy based on the C e g method of Ref. [14| . The one-cx uncertainty and the 
0.3keVee analysis threshold are also shown. 
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Figure 2: Energy scale comparison between calibration data and simulation results. 
Calibrations were conducted at 122keV ( 57 Co), 88keV, 58keV, 22keV, 8keV ( 109 Cd), 
59.5keV, 17.8 keV ( 241 Am), and 5.9keV ( 55 Fe). The error bars in the figure stem from 
uncertainties in the surface reflection properties of the calibration source housing within 
the simulation. 
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Figure 3: Expected trigger efficiencies for 7GeV, 12GeV, and 18GeV WIMPs derived 
from LED data in the XMASS detector. 
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Figure 4: The Cherenkov event (head-to-total ratio > 0.6) cut for the XMASS data (left) 
and for simulated WIMP events (right) as a function of the observed energy. 
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Figure 5: The obtained energy spectrum at each step of the data reduction. Raw data is 
shown as the solid line, the result of the 10 ms to the previous event timing cut appears as 
the dashed line, and the dash-dotted line shows the data after application of the timing 
distribution width cut. The filled points show the result of the Cherenkov event cut. 
Details are presented in the text. 
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Figure 7: Simulated WIMP energy spectra in the XMASS detector assuming the maximum 
cross section that provides a signal rate no larger than the observation in any bin above 
0.3 keVee. 
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Figure 8: Spin- independent elastic WIMP-nucleon cross section as a function of WIMP 
mass. All systematic uncertainties except that from £ e ff are taken into account in the 
XMASS 90% C.L. limit line. The effect of the C c g uncertainty on the limit is shown in 
the band. Limits from other experiments and favored regions are also shown [3, [US SHE]- 
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